Majorana zero modes in graphene with trigonal warping 
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We study the low energy properties of warped monolayer graphene, where the symmetry of the 
original honeycomb lattice reveals itself. The zero energy solutions are Majorana fermions, whose 
wavefunction, originating from the corresponding modified Dirac equation is spatially localized. 
Experimental consequences are discussed. 
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Graphene, a single sheet of carbon atoms have at- 
tracted enormous attention since its discovery in 20041]. 
Due to the interaction of electrons with the background 
C-atoms, its quasiparticles obey to the Dirac equation in- 
stead of the Schrodinger one, with the velocity of light (c) 
replaced by the Fermi velocity (vp ~ c/300) in graphene. 
This gives rise to many interesting physical phenom- 
ena, including the unconventional Hall effect, Klein tun- 
neling, Zitterbewegung, and peculiar electric and heat 
conduction 

In realistic graphene, there are several physical effects 
that appear, which can have a significant impact on its 
properties. The most important of these effects is the in- 
fluence of ripples [§], leading to weak randomization, as 
scattering of short-range defects do not conserve isospin 
and valley [4j and trigonal warping [5J of the electronic 
band structure which introduces asymmetry in the shape 
of the Fermi surface about each valley, as have been mea- 
sured experimentally [(|. 

These two effects will also strongly alter the chirality 
of the charge carriers in graphene [|, 0|. Ripples, i.e., 
disorder leads, to a description of the system in terms of 
a Grassmannian Non Linear Sigma Model belonging to 
the symplectic symmetry class [7| ■ Trigonal Warping, on 
the other hand, breaks [B| the Time reversal Symmetry 
[j| and thus breaks the supersymmetry of the Non Linear 
Sigma Model, leading to the appearance of Majorana zero 
modes. 

While much attention has been channeled into the 
study of ripples, there is only a few approaches to realis- 
tic trigonal warping effects. Hence, in this letter we want 
to fill this gap, by studying Majorana zero modes due to 
trigonal warping. An early computation of the effective 
action as sketched in 0] shows that - with trigonal warp- 
ing - the Non Linear Sigma Model belongs to the unitary 
symmetry class. This induces [§] the disappearance of 
the perfectly conducting channel in a way analogous to 
what happens when a magnetic field is applied. 

Hereafter, contrary to Ref. 0], we derive the Majo- 
rana zero modes from the Dirac equation, by writing the 
Hamiltonian, e.g., in K point, by explicitly adding a trig- 



onal warping term j3(z)d1 
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iky with the well-known 



Here z = k 

notations k x — —i(d/dx) and k y = —i(d/dy). 

The trigonal warping term (3(z) = (3(r) exp(3i£8) is a 
generalized form of the second order contribution to the 
spectrum derived in Ref. [5|: 
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where a is the lattice constant, j3 is an adjustable param- 
eter, of the order of unity [H, for which we introduce a 
slow r variation, i.e., (3 = 0{r), however we will see that 
the final result does not depend on this choice, but rather 
by the vorticity given by exjp(3i£6). 

The new Hamiltonian term from Eq. (J2]) is not in- 
variant under the time reversal operation. However, it 
should be noted that the total Hamiltonian given in Eq. 
([T]) does not break the original time reversal symmetry. 
In fact K and K' are related to each other by time re- 
versal. If there is not intervalley scattering between K 
and K', the original time-reversal symmetry is irrelevant 
anyhow. 

The wave-function in the K point we write as 



A" 



(3) 



hence the uniform limit gives back the well-known energy 
spectrum: 



£±{kx, ky) ~ {fc 2 H 
[(fc 2 - 3kl)k x cos3(9+(fc 2 



k y 



2,k 2 y )ky sinW]} 1 / 2 



(4) 



As it can be seen, the second order term gives rise, 
indeed, to a trigonal warping of the dispersion, where 
the upper and lower sign corresponds to the conduction 
and valence bands, respectively. 
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FIG. 1: (Color online) The first Brillouin-zone of the honey- 
comb lattice (dashed lines) with the orientations of the co- 
ordinate axis and the equienergy surfaces close to the Dirac 
point with trigonal warping are shown. 



Going beyond the uniform limit, the full solution of 
the Dirac equation takes the form: 



—id z u = Eu , 
+id?v = Ev . 



(5) 



We observe that for £ even, I — 2n , we can completely 
eliminate the phase dependence from Eq. ([5]) by making 
a transformation: 



ue 



3/ n 



ve 



-3in9 



(6) 



Hence, Eq. are topologically equivalent to the uni- 
form limit, presented just earlier. 

The situation is completely different if I is odd, t = 
2n — 1. In this case we cannot eliminate the phase de- 
pendence. Even with the help of a transformation ^ a 
phase factor exp(3i6>/2) always remains in Eq. due 
to which, it can be easily verified that, e.g., the complex 
conjugate of the second equation from ([5]), i.e., is iden- 
tical to the first one, u. This guarantees one localized zero 
mode with u — u*. 

In the next step, we attempt to obtain the zero modes 
explicitly. For this we write the Dirac equations ^ in 



compact form as: 

- id z u(r, 9) + VW)d- z (VP&)d- z ) u(r, 6) = . (7) 
Which, in polar coordinates becomes: 
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u(r,6) = 0, (8) 



where, as mentioned earlier, we have odd vorticity, I = 
In — 1, and for simplicity we used the notation 

F(r,6) = f(r)e me / 2 e- ie , (9) 

with (3(r) = f 2 (r). The radial solution of Eq.® is similar 
to the zero modes studied by Ghaemi and Wilczek [101 ] : 



fu"+ 1 + ff 



1/2 a 



, (10) 



The solution of which for large distances, r 3> a, 
ponentially decreasing 

u w exp(-r// 2 ) . 

While for r«a tends to a constant value 



exp 



1/2 / ; 



(11) 



(12) 



hence independent of the form of f(r), i.e., (3{r). 

A similar solution can be obtained for K' hence, we 
have two localized Majorana zero modes in both K and 
K'. The number of solutions for K equations minus the 
number of solutions for K' equations equal the vorticity. 
This is the index theorem 11( already stated f or g raphene 
with geometrically induced vorticity in Ref. |12j |. 

This is a necessary requirement in order to have un- 
paired, real Majorana modes. In spite of the beautiful 
simplicity in deriving Majorana fcrmions, they are not 
easy to come by in nature: if the pairs of Majorana modes 
are not localized and kept apart they will immediately 
form a Dirac fermion. And more stringently, any system 
to exhibit Majorana modes has to be a perfect diamag- 
net, because the interaction of the complex fermions with 
photons is not diagonal. The real and imaginary compo- 
nents (which are the Majorana fermions) would rapidly 
mix by electromagnetic interaction. Graphene is know 
to be a diamagnet [Hj], hence Majorana fermions could 
exist. 

However, other systems where Majorana could be seen 
has to comply to these requirements. The easiest way 
to construct a perfect diamagnet is through a supercon- 
ducting state, and hence other systems where Majorana 
zero modes have been shown to exist are superconduct- 
ing states, such as the chiral spinless two-dimensional su- 
perconductors characterized by p x + ip y order parameter 
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14 1 and topological insulator/superconductor structures 
151 ]: or fractional quantum Hall states, such as the the 
v — 5/2 state One of the main difference between 

these systems in view of the Majorana modes is that in 
the v = 5/2 quantum Hall state, a Majorana bound state 
is associated with charge e/4 quasiparticle. While in all 
other system, including graphene, the Majorana modes 
are electrically neutral. In the fractional quantum Hall 
state, thanks to the e/4 charge the quasiparticle's non- 
Abelian statistics can be probed by measuring charge 
transport of the edge states 17J. And even more di- 
rectly recent experiments [l8j have shown evidence for 
the quasiparticle charge. However, detecting the Ma- 
jorana fermions in superconductors or graphene is still 
lacking direct experimental confirmation. 

Indirect confirmation of Majorana fermions in 
graphene may emerge from already well-known experi- 
mental findings. In this respect it is interesting to re- 
mark, firstly, that zero modes resulting from ripples, i.e., 
bond, charge, i.e., potential disorder [191 ] or random vec- 
tor potential, have different properties from our solu- 
tion, namely disorder of these types will suppress electron 
backscattering, this is known as the anti-localization ef- 
fect [20I ]. In our case, however, there is no suppression of 
backscattering, which leads to conventional weak local- 
ization. Hence, if weak localization is seen to be domi- 
nating in graphene, this may well be due to topological 
terms of Majorana type. 

The second possibility is related to the minimal con- 
ductivity of graphene [2l| . It is believed that an intuitive 
explanation of this phenomenon is that, since chirality is 
due to the equivalence of the K and K' and produces 
nonequivalent Dirac cones at both corners of the Bril- 
louin zone, the wave functions of electrons have an isospin 
quantum label and, thus, in isospin conserving scattering 
process, i.e., processes unable to distinguish between K 
and K', charge carriers cannot be backscattered and the 
resistance is reduced. 

When the Non Linear Sigma Model is used to de- 
scribe the transport of two-dimensional Dirac electrons 
in a random electrostatic potential, a topological term is 
always required in the mathematical formulation 0, 1221 ] . 
The topological term arises from Majorana fermions, and 
hence the presence of Majorana zero modes can change 
low energy (long distance) properties of the graphene 
drastically and can be responsible for its transport prop- 
erties. 

In conclusions, we have shown that taking into account 
trigonal warping explicitly, we can account for the pres- 
ence of Majorana fermions in single layer graphene. Ex- 



perimental consequences are discussed. 
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